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Abstract 

Hydrodynamic turbulence is studied as a constrained system from 
the point of view of metafluid dynamics. We present a Lagrangian 
description for this new theory of turbulence inspired from the analogy 
with electromagnetism. Consequently it is a gauge theory. This new 
approach to the study of turbulence tends to renew the optimism to 
solve the difficult problem of turbulence. As a constrained system, 
turbulence is studied in the Dirac and Faddeev-Jackiw formalisms 
giving the Dirac brackets. An important result is that we show these 
brackets are the same in and out of the inertial range, giving the way 
to quantize turbulence. 
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1 Introduction 



Since its discovery, turbulence has attracted an enormous effort from sci- 
entists to study it H, |]. In recent years, there has been an increasing 
interest in studying turbulent hydrodynamics because lots of phenomena in 
Nature are turbulent, like in the flow of blood in biomechanics, meteorology 
and ocean engineering, astrophysics and formation of galaxies in the Uni- 
verse. Beyond these phenomena, turbulence has became a very fruitful field 
of research for theoreticians, who have studied the analogies between tur- 
bulence and field theory, critical phenomena and condensed matter physics 
[^, 1^, ^, ^ ^, ^, |TT], leading to renewed optimism to solve the problem 
of turbulence described by the Navier-Stokes equations 

+ vW^u, (1) 

where u{x,t) is the velocity field, w{x,t) the vorticity field, p{x,t) is the 
pressure, p is the density and u the kinematic viscosity. 

Among turbulent phenomena, the hydrodynamic turbulence is one in 
which physicists have enormous interest because of the universal character- 
istics stressed by an incompressible fluid at high Reynolds numbers in the 
fully-developed turbulent regime. The regime related to the limit of high 
Reynolds number, R ^ oo, with R = {LV)/^, which measures the competi- 
tion between convective and diffusive processes in an incompressible fluid de- 
scribed by the Navier-Stokes equations. Above, L is the integral length-scale 
of the largest eddies that should appear and V, is a characteristic large-scale 
velocity. 

An important feature in turbulence is the presence of vorticities , but 




another flow quantity, the Lamb vector, can be used to describe turbulence 114 



together with vorticities W5\ in the so called metafluid dynamics. In metafiuid 
dynamics, the study of average quantities of an incompressible fluid in the 
fully developed regime is proposed, given a system where the average fields 
show up in a continuum inter-relation and respond as waves to quantities 
named turbulent sources. To do so, the Lamb vector / and the vorticity be- 
come the kernel of the turbulent dynamics, being related by l{x,t) = w x u. 
Defining Bernoulli energy function 0, (f){x,t) = - + one can introduce a 
theoretical concept, the turbulent charge density n, n{x,t) = — V^0. In the 
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inertial range, the Bernoulli energy function is conserved, the energy gener- 
ated at larger scales is transferred to smaller wave numbers across the region 
and energy pumping and dissipation are not relevant allowing us to make a 
very close analogy between turbulent hydrodynamics and electromagnetism 
[p^ . Later on, we show there is no loss of generality in treating the turbulent 
flow in inertial range from the point of view of our proposal. 

This analogy between hydrodynamic turbulence and the electromagnetism 
in inertial range makes possible to study metafluid dynamics. So as in 
the electromagnetism |]l6l, where we write the Lagrangian density as £ = 
^{E'^ — B'^) {E is the electric field and B is the magnetic field), one can write 
down the Lagrangian density for the theory of turbulence as £ = |(/^ — cJ^) 
or £ = -lF^^F^"'-J^V'', with F^"' = d'^V - d^V^", = (0, m), and J^' = 
{n, J), J = un + V X {u.uj)u + uj x V(0 + m^) + 2(/.V)m is the turbulent 
current density. Writting J = Jtr + Ji, where Jtr is the turbulent transverse 
current density and J; = V0 is the turbulent longitudinal current density, 
the Lagrangian density becomes C = -iF^^F^^ - J^V^' + where now 

= (n, Ji). We note that V. Jtr = 0. 

Despite of all resemblance between hydrodynamic turbulence and the 
dynamics of electromagnetic fields, there is a conceptual difference in the 
identification of the physical entities. In the classical electromagnetism, the 
physical fields are the electric and magnetic fields and the potentials are just 
mathematical artifices while in the metafluid dynamics the potentials play 
physical role. 

The analogy between turbulent hydrodynamics and electromagnetism 
suggests we can study the turbulence as a constrained system. The the- 
ory of constrained systems plays an important role in the study of physical 
systems. It includes almost all known fundamental theories. Constrained 
systems were first studied systematically by Dirac p7[] . 

Constrained systems are characterised in phase space by the presence of 
constraints, which are functions of the coordinates and momenta. These 
constraints can be classified into primary, secondary, etc., or first and second 
class [O, [TB|, ITPI]. First class constraints imply the presence of gauge invari- 



ance in the theory, since they generate gauge transformations. Beside that, 
they exhibit the structure of the corresponding gauge group. On the other 
hand, second class constraints do not have these properties. 

There exist several methods to treat constrained systems based on the 
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classification given above. Most of tliese deal with first class constraints, for 
the second class ones there is the Dirac bracket method. There is also a more 
recent method by Faddeev and Jackiw [^], which does not follow the above 
classification. In this method, the Dirac formalism can be avoided and basic 
bracket relations can be obtained without using the usual Dirac brackets. 

We will consider the Dirac and Faddeev- Jackiw formalisms to treat tur- 
bulence as constrained system. Those procedures have been used with great 



success in Quantum Field Theory for quantizing some models[21, 22 



2 Dirac formalism 



Under the point of view of a geometric interpretation, Arnold [^] showed. 



using Lie algebra, that Euler flow can be described in the Hamiltonian for- 
malism in any dimension. This has a lot of interesting consequencies for fluid 



mechanics and has been studied ^ , but it is not quite obvious that this 



method could be used for viscous fluid. That is where the metafiuid dynam- 
ics comes, showing a way to find a Hamiltonian formalism even for turbulent 
flow with viscosity. 

Using again the analogy between electromagnetism and turbulence, one 
obtains the coniugate momentum for £ as tt'^ = 4^ = that immedi- 

ately leads us to a primary constraint ^ 0. So^ the primary Hamiltonian 
becomes Hp = J dx {^T^'V^ - £ + AttO) = J dx (ivf^ + i(V x m)^ + 27r.V0 

In order to apply the Dirac Hamiltonian method [|T7|, we need to look for 
secondary constraints. Imposing that primary constraint must be conserved 
in time, we get = {vr^. Hp} ^ 0. This consistency condition over the 
primary constraint leads to the secondary constraint V.vr — J° ^ 0, that is 
exactly one of the metafiuid dynamics equations. If one goes on and imposes 
the consistency condition over this secondary constraint, we observe that no 
new constraints are obtained via this iterative procedure. 

The constraints 7r° «i and V.vf — J° ^ are first class constraints. It 
means this new theory of turbulence is gauge- invariant. A result expected 
due to the analogy between electromagnetism and turbulence. This anal- 
ogy allows us to apply to this new theory of turbulence all the machinery 
associated with gauge theories. Since we have two constraints of first class. 
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we have two degrees of freedom to field V^. These two degrees of free- 
dom can be associated to the vorticity. This means we must explore the 
gauge invariance of the theory to choose two gauges. Following the anal- 
ogy between electromagnetism and turbulence we find the gauges V.m = 
0, (Coulomb or transverse gauge) and = a (constant) (Temporal gauge), 
where the first gauge comes from condition of incompressibility of fluid and 
the second, from the fact that Bernoulli energy function is constant. 

So, the theory has the constraints Ti = — a ^ 0, T2 = tt^ ^ 0, T3 = 
V.M ~ 0, and T4 = V.vf — ^ 0,which are second class ones. To implement 
the Dirac brackets we need to compute the matrix elements of their Poisson 
brackets, which reads Cu = S{x — y) = — C21, C34 = —V^6{x — y) = 
— C43, and all other elements are zero. Using the Dirac brackets {A, B}* = 
{A, B} - {A, T^}C-}{Tf,, B}, we get 

{ui{x),Uj{y)}* 
{ui{x),'Kj{y)}* 

Inspite of these Dirac brackets have been computed in the inertial range, 
they are valid for all fully-developed turbulence, as we show bellow. 



= = {iTiix), TTjiy)}*, 

= U,-^U{S~y). 



3 Faddeev-Jackiw formalism 

With the purpose to show that the presence of viscous term in the Navier- 
Stokes equation does not change the Dirac brackets, let's analyze the metafiuid 
dynamics using a second treatment for this constrained theory, the Faddeev- 
Jackiw procedure, also known as sympletic formalism. 

Faddeev and Jackiw |^ showed how to implement the constraints di- 
rectly into the canonical part of the first order Lagrangian. Once this method 
is applied to first order Lagrangians, for implementing the Faddeev-Jackiw 
procedure we write the Lagrangian density of turbulence making explicit the 
time-derivative of the fields : 

1 • • 1 1 

C = _^2 + M.V0+ -(V0)^ - -(V X m)2 - z/u.V^u 

- u(y^u).{V(p) + ^u\V^uf + U.J- uu.Vn- (pn, (3) 
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and for transforming it from second to first order, we consider the momentum 
as the auxihary field, extending the configuration space. In that way, we get 
the equation of motion for vr : vf = — V0 — u + z/V^m. Replacing back this 
result in the Lagrangian we return to the original quadratic term. Using that 
J = Jtr + Jh the Lagrangian density in eq.(|^) becomes 

C = -iyf^ - 7f.V0 - vf.M + z/TT.V^n + ^(V X u)^ 

— uu.Vn + u.Jtr + u-Vc/) — ncj) (4) 

In this case, the first order Lagrangian reads C^^^ = —n.u — u. 
where the potential density is f/^^^ = |7f^ + 7f.V0 — |(V x ■u)^ — ut^.V'^u — 
u.JtT + (t)n + uu.Vn. 

The initial set of sympletic variables ^f''^ = {ui, iTi, 0} allows us to 
identify the quantities a-°^" = — vTj — af^'^ = 0, and af^^ = 0. These 
lead us to the following matrix f^^\ 



/W= df 6,, \s{x-y), (5) 




which obviously is singular. Thus, the system has constraints in the Faddeev- 
Jackiw formalism. The components of the zero-mode v^^^ = (f*^, 0, v'^) sat- 
isfy vf - diV^ = 0. 

The primary constraint is obtained from / dxf "^(x) [— V.tt (x) + n{x)] = 0. 
Since v'^{x) is an arbitrary function and the fluid is taken incompressible, we 
obtain the constraint V.7f(x) — n{x) = 0. 

Introducing this constraint back into the Lagrangian by means of a La- 
grange multiplier A one obtains C^^^ = —n.u — V0.u -|- A(V.7f — n) 
where U^^^ = ^n^ — |(V x — z/tt.V^m — u.Jtr + yu.Vn. 

Considering now that the new set of sympletic variables is given in the 
following order, ^ = (-Uj, 0, vTj, A), we have af^'^ = —iTi — di(j), af'^'^ = 
0, and a*^^-''** = V.vf — n. From there we obtain the sympletic matrix f^^\ 

/ Of 6,, \ 
V -6»J / 
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that is a nonsingular matrix. Then we can identify it as the sympletic tensor 
of the constrained theory. The inverse of Z^^-' gives us the Dirac brackets of 
the physical fields and can be obtained in a straightforward calculation. The 
result is 



(1)^-1 





_A. . J ' J 










1 

V2 



A. . ' J 



From the matrix (Z^^-*) ^ we identify 





1 

V2 
V2 





5{x-y) 



(7) 



(8) 



where the auxiliary fields were eliminated and we returned to the original 
variables. These brackets are nothing more than brackets of the turbulence 
gauge field in the gauge V.m = 0, obtained in eq. (|^), showing us that the 
Dirac brackets are the same for the Navier-Stokes equation with or without 
the viscous term. 

The straightforward step would be the quantization of this constrained 
theory, that can be obtained by the well known canonical quantization rule 
{ , }* ^ —i[ , ]■ Doing so, we get the commutators 



[Ui[X),Uj 



[Ui{X),TTj 




6{x-y). 



(9) 



Once we have the canonical quantization rule, we can apply standard 
Quantum Field Theory technics to find the generating functional, conse- 
quently correlation functions and all physical quantities [^, one wish. 



4 Conclusion 



In conclusion, we have proposed to study incompressible turbulent hydrody- 
namics in the metafiuid dynamics formalism applying tools commonly used 
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in Quantum Field Theory for constrained systems. This study is possible be- 
cause of the close analogy between incompressible turbulent hydrodynamics 
and the electromagnetism. Exploiting this analogy we computed the Dirac 
brackets using the Dirac and Fadeev-Jackiw formalisms. We showed that 
Dirac brackets found in the Dirac formalism in the inertial regime {u = 0) is 
valid for all regime of fully-developed turbulence, as we found in the Fadeev- 
Jackiw formalism. We believe this approach of turbulence as a gauge theory 
can lead to results that tend to renew the optimism to solve the difficult 
problem of turbulence. 
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